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We investigate the source of large variations in the observed detection efficiencies of supercon- 
ducting nanowire single-photon detectors between many nominally identical devices. Through both 
electrical and optical measurements, we infer that these variations arise from "constrictions:" highly 
localized regions of the nanowires where the effective cross-sectional area for superconducting cur- 
rent is reduced. These constrictions limit the bias current density to well below its critical value 
over the remainder of the wire, and thus prevent the detection efficiency from reaching the high 
values that occur in these devices only when they are biased near the critical current density. 



PACS numbers: 74.76.Db, 85.25.-j 



Superconducting nanowire single-photon detectors 
(SNSPDs) (HUSH provide a unique combination of 
high infrared detection efficiency (up to 57% at 1550nm 
has been demonstrated) and high speed (<30 ps tim- 
ing resolution 0, [B[ , and few-ns reset times after a detec- 
tion event Q). Applications for these devices already be- 
ing pursued include high data-rate interplanetary optical 
communications [f|, spectroscopy of ultrafast quantum 
phenomena in biological and solid-state physics 0, @|, 
quantum key distribution (QKD ) [9[ , and noninvasive, 
high-speed digital circuit testing flO ]. 

In many of these applications, large arrays of SNSPDs 
would be extremely important []|. For example, exist- 
ing SNSPDs have very small active areas, ma king optical 
coupling relatively difficult and inefficient @, hlfl- Their 
small size also limits the number of optical modes they 
can collect, which is critical in free-space applications 
where photons are distributed over many modes, such 
as laser communication through the atmosphere (where 
turbulence distorts the optical wavefront) and in fluores- 
cence detection. Furthermore, it was shown in Ref. [4[ 
that the maximum count rate for an individual SNSPD 
decreases as its active area is increased, due to its kinetic 
inductance, forcing a tradeoff between active area and 
high count rates. Count rate limitations are particularly 
important in optical communications and QKD, affect- 
ing the achievable receiver sensitivity or data rate [f| fl2| . 
Detector arrays could provide a solution to these prob- 
lems, giving larger active areas while simultaneously in- 
creasing the maximum count rate by distributing the flux 
over many smaller (and therefore faster) pixels. Large ar- 
rays could also provide spatial and photon-number reso- 
lution. Although few-pixel detectors have been demon- 
strated [1, [i| [lj, [lH, fabrication and readout methods 
scalable to large arrays have not yet been discussed. 

A first step towards producing large arrays of SNSPDs 



is to understand (and reduce) the large observed varia- 
tion of detection efficiencies for nominally identical de- 
vices [2(|, which would set a crippling limit on the yield 
of efficient arrays of any technologically interesting size. 
In this Letter, we demonstrate that these detection effi- 
ciency variations can be understood in terms of what we 
call "constrictions:" highly localized, essentially pointlike 
regions where the nanowire cross-section is effectively re- 
duced, and which are not due to lithographic patterning 
errors (line-edge roughness). 

The electrical operation of these detectors has been 
discussed previously by several authors [J, 0, d, 0, [T^ . 
so we only summarize it here. The NbN nanowires are 
biased with a DC current /bias slightly below the crit- 
ical value Iq- An incident photon of sufficient energy 
can produce a resistive "hotspot" which in turn disrupts 
the superconductivity across the wire, producing a series 
resistance which then expands in size due to Joule heat- 
ing 0, EH- The series resistance quickly becomes ^>50il, 
and the current is diverted out of the device and into the 
50J1 transmission line connected across it, resulting in a 
propagating voltage pulse on the line. The device can 
then cool back down into the superconducting state, and 
the current through it recovers with the time constant 
Lk/50f2, where L\ is the kinetic inductance 

The nanowires used in this work were patterned at 
the shared scanning-electron-beam-lithography facility in 
the MIT Research Laborator y of Electronics using a pro- 
cess described in Refs. 0,13:13' on ultrathin (~ 5 nm) 
NbN films grown at Moscow State Pedagogical Univer- 
sity [15| . The majority of the wires were on average 90 
nm in width, and were fabricated in a meander pattern 
with 200 urn line pitch (45% fill factor), subtending an 
active area of 3 x 3.3 urn or 10 x 10 urn. Some devices 
with 54 nm average width and 150 nm pitch (36% fill fac- 
tor) were also fabricated - see Fig. . The devices had 
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FIG. 1: Figure [T] (color online) Variations in SNSPD detec- 
tion efficiency, (a) Histogram of the DEs measured for 132 
devices from a single fabrication run, on a single chip. The 
devices were 3 x 3.3 [im in size, and composed of a 50 urn 
long nanowire in a meander pattern with 45% fill factor. The 
measurements were made at T — 1.8K, and /bias = 0.975/c bs . 
Note that the peak DE of 22% increases to 57% with the 
addition of an optical cavity, as described in Ref. |2J (b) 
Measurements of the DE vs. Ibi^a/Ic* 3 for a selection of these 
devices. 



critical temperatures Tq ~ 9 — 10 K, and critical current 
densities J c ~5x 10 10 A/m 2 at T = 1.8K. 

The experiments were performed at MIT Lincoln Lab- 
oratory, using the procedures and apparatus discussed in 
detail in Refs. 0, H[. Briefly, the devices were cooled to 
as low as 1.8 K inside a cryogenic probing station. Elec- 
trical contact was established using a cooled 50 mi- 
crowave probe attached to a micromanipulator, and con- 
nected via coaxial cable to the room-temperature elec- 
tronics. We counted electrical pulses from the detectors 
using low-noise amplifiers and a gated pulse counter. To 
optically probe the devices, we used a 1550 nm mode- 
locked fiber laser (with a 10 MHz pulse repetition rate 
and <1 ps pulse duration) that was attenuated and sent 
into the probing station via an optical fiber. The devices 
were illuminated from the back (through the sapphire 
substrate) using a lens attached to the end of the fiber 
which was mounted to an automated micromanipulator. 
The focal spot had a measured e -2 radius of ~ 25 u.m. 

Figure Q] illustrates the detection efficiency (DE) vari- 
ations observed on a single chip of 132 devices of the 
same geometry, fabricated in a single run. These devices 
are the same ones reported in Ref. Q before the optical 
cavities were added (maximum DE after the addition of 
cavities was 57%). In panel (a) we show a histogram of 
the measured detection efficiencies at /bias = 0.975/g bs 
(where ig bs is the observed critical current of each de- 
vice), and (b) shows some representative data of the ob- 
served DE as a function of /bias/-fc bs - Note that the shape 
of these curves also varies significantly. As we show be- 
low, these data can be explained with the hypothesis, 
first suggested in Ref. , that some devices have "con- 
strictions:" regions where the (superconducting) cross- 
sectional area A cs of the wire is reduced by a factor we 
label C . This effectively reduces the observed critical 
current by that same factor (Ig bs = Jc^csC = IcC), 



and prevents the current density everywhere but near the 
constriction from ever approaching the critical value (and 
hence prevents the wire from having a high DE except 
locally near the constriction). 

If all the nanowires were identical in all dimensions save 
for a pointlike constriction, we would expect that if the 
data of Fig. [ljb) were plotted vs. absolute current /bias 
(rather than /bias/-fc bs ) it would all lie on a single, uni- 
versal curve, with the data for more constricted devices 
simply not extending to as high currents. This turns out 
to be approximately true, but variations from device to 
device across the chip either in film thickness or in the 
nanowire width obscure this essential feature of the data 
when it is plotted in this simple way. Instead, we present 
our results as shown in Fig. [2] In panel (a) , the DE data 
for each of 170 devices with 90 nm wide wires (across two 
chips fabricated in separate runs) is shown superposed 
(filled circles indicate T = 1.8K, and crosses T — 4.2K). 
All of the data for each temperature can be made to lie on 
a single universal curve by scaling the critical current I£ bs 
for each device by an adjustable factor (which is just 1/C: 
Ic = ^c bs /C)- The very fact that data from this many 
devices can be so well superposed is already suggestive of 
a universal shape. However, we can now take the C value 
for each device extracted using this scaling procedure and 
cross-check it. Based on our previous discussion, if all 
wires were identical save for constrictions, we would ex- 
pect C = I^ s /Ic, i.e. C to be exactly proportional to 
7g bs . Due to the variations across a chip discussed above, 
this is only partially true. However, we can normalize out 
these variations using a very simple method: instead of 
comparing C directly to I£ bs , we instead compare it to 
the product Ig bs i? n = (C x J c A cs )(pJ/A cs ) = C x J c pJ, 
where A cs and R n are the cross-sectional area and room- 
temperature resistance of each nanowire, Jq and p n are 
the critical current density and room-temperature resis- 
tivity of NbN, respectively, and I is the total wire length. 
This product depends on the wire geometry only through 
I (which is fixed lithographically and does not vary appre- 
ciably between devices) and not on each wire's individual 
A cs . Figures HJb) and (c) show a comparison between the 
C values extracted from the data in Fig. [Ha) and the 
7g bs i? n product. The data lie on a straight line through 
the origin, indicating that these two independent mea- 
sures of C are mutually self-consistent. 

Our data can also be used to infer that the constric- 
tions are essentially pointlike (i.e. very short in length). 
The open circles in Fig. HJa) are data for 15 devices with 
54 nm wide wires, and clearly show a dramatically differ- 
ent shape (i.e. high DE persists to much lower currents 
than for the wider wires). The broken lines are estimates, 
based on the data for 90 nm and 54 nm wide wires, of 
what one would expect for a device having 90 nm wide 
wire, except at a single constriction 54 nm wide (cor- 
responding to C ~ 0.6) with a length of either 0.5 urn 
(dashed line) or 50 nm (dotted line) long. These curves 
have a different shape from the data for 90 nm wide wires, 
because at low currents the DE is dominated by the re- 
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FIG. 2: Figure[2] (color online) Constriction values extracted 
using DE vs. /bias data, (a) Universal DE curves for devices 
with 90nm wide wires, at T = 1.8K (•) and T = 4.2K (+). 
Data from 170 devices distributed over two separate fabrica- 
tion runs is shown superposed. By rescaling the Iq 13 of each 
device such that all data lies on a single universal curve as 
shown, the constriction C (which indicates the fraction of the 
wire's cross-sectional area that is superconducting at the con- 
stricted point) can be obtained. Also shown (o) are data for 
15 devices having 54 nm wide wires (and 36% fill factor) at 
T — 1.8K, indicating that narrower wires exhibit a very dif- 
ferent universal curve shape. These data provide evidence for 
the localized nature of the constrictions, since any appreciably 
long section of wire having a smaller cross-section should sig- 
nificantly alter the shape of the curves, making it impossible 
to superpose them as shown. This is illustrated by the broken 
lines, which show a simple estimate for a 90 nm wide wire hav- 
ing a 54 nm wide constriction with a length of 0.5 urn (dashed 
line) or 50 nm (dotted line). These estimates are obtained 
simply by adding together the universal curves for the two 
wire widths, in a ratio given by the length of wire with each 
width, i.e. (0.5 (J.m/50 nm) and (0.05 [J.m/50 urn). Note that 
these broken lines terminate at Ibiaa/Ic" 1 x C = 0.6 = 54/90 
due to the assumed 54 nm constriction of the 90 nm-wide 
wire, (b) and (c) Constriction values C obtained from the 
data in (a) (for (b) T = 1.8K and (c) T = 4.2K) vs. those 
obtained using the IqRxx product. The C values in both 
cases are normalized absolutely using Lk^bias), as described 
below (see Fig. The solid lines are straight lines through 
the origin with slope 1; no fitting was used. 

gion of wire near the constriction, while at higher current 
it becomes dominated by the contribution from the rest 
of the wire length. This very different shape should be 
distinguishable if it were present, and should prevent the 
data from being superposed onto a single curve. The ab- 
sence of this in our data at any level above the noise indi- 
cates that the constricted regions are likely much shorter 
than ~0.5u-m. 

So far, we have in fact only measured constriction in 
a relative sense; that is, we have no way to tell if our 
best devices in fact have C = 1. To address this, we can 
exploit the known dependence of the kinetic inductance 
on bias current. Kinetic inductance arises from energy 
stored in the effectively ballistic motion of Cooper pairs; 
as the current density is increased towards the critical 
value, the density of Cooper pairs is depleted, forcing 
the remaining pairs to speed up (and therefore store more 
kinetic energy per unit volume) to maintain the current. 



Hence, the kinetic inductance increases as J — * Jc [17j . 
Since the kinetic inductivity locally increases with J/ Jc, 
the total kinetic inductance of the wire provides a way to 
determine if the current density is indeed near the crit- 
ical value over the whole wire or only at one localized 
place. We measured the inductance of our nanowires 
using a network analyzer, by observing the phase of a re- 
flected microwave signal as a function of frequency. We 
then fit this data using a suitable electrical model to ex- 
tract the inductance value. A bias tee was inserted into 
the signal path to superpose the desired /bias with the 
network analyzer output. The phase contributions from 
the coaxial cable, bias tee, and microwave probe were re- 
moved by probing an in situ microwave calibration stan- 
dard (GGB Industries CS-5). The microwave power used 
in this measurement corresponds to a peak current am- 
plitude of < 0.5u^4, and the critical current measured in 
the presence of the microwaves was within 10% of that 
measured in their absence (~ 20|jA for typical devices at 
T=1.8K). 

In Fig. 03a), we show the measured inductance vs. 
current of two devices; one which has nearly the highest 
detection efficiency observed on this chip (22% - filled 
circles), and the other having one of the lowest (0.1% - 
filled triangles). Also shown is the prediction for I/k(^bias) 
from Ginsburg-Landau theory (see, e.g., |17|). with no 
free parameters (solid line). The data for the high-DE 
device show good agreement with this prediction, indi- 
cating that this device is indeed unconstricted. However, 
for the low-DE device the inductance does not increase 
as much as predicted, which is precisely what we would 
expect within the constriction hypothesis; the current 
density is only near critical at one localized place (which 
constitutes a negligible fraction of the total wire length) 
whereas everywhere else the current density is lower, pro- 
ducing a smaller increase in inductance. The factor by 
which Jg bs must be rescaled for a given device so that 
the ik(^bias) data matches the prediction constitutes an 
absolute measurement of C. One such measurement for 
a single representative device, from among a large set of 
nominally identical devices, then allows us to correctly 
normalize the C values obtained using either of the pre- 
vious methods described above (see Fig. [2]) for all other 
devices in that set. 

We can also verify that the observed ik(^bias) and 
I c hs R n product give mutually consistent results. To 
check this, for each device we measured the inductance 
ratio TZ L = L(I o )/L(0), where J « 0.9/g bs . Using the C 
obtained from the normalized /g bs i? n product, we also 
obtain Iq/Iq — IqC/I^ 5 for each device. The inset 
to Fig. [3|a) shows TZl vs. Io/Ic (filled squares), and 
the data are in reasonable agreement with the Ginsburg- 
Landau prediction (solid line). 

In addition to providing evidence for the constriction 
hypothesis, the measurement of Lk(^bias) and Jg bs i? n pro- 
vide a powerful diagnostic tool, since they constitute a 
purely electrical measurement of C, which can then be 
used to predict the detection efficiency, as indicated in 
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FIG. 3: Figure [3] (color online) Absolute measurement of 
constriction using bias current-dependence of kinetic induc- 
tance, (a) The kinetic inductance of our nanowires should 
increase with /bias, due to the depletion of the Cooper-pair 
density near the critical current density (solid line). A de- 
tector with the highest observed detection efficiency on the 
chip (22%) behaves as expected (•), with no free parameters. 
However, a detector with much lower DE (0.1%) does not 
(filled triangles). This is due to a constriction, which pre- 
vents the current density from approaching critical anywhere 
but near this one localized point. Inset: the inductance ratio 
TZl = L(I » 0.9Ig bs )/L(0) measured for each device (filled 
squares), plotted vs. Iq/^I^/C) = Io/Ic, where C is ob- 
tained from the I^Rn product. These data agree with the 
prediction (solid line), indicating that IZl and I^Rn give 
mutually consistent results for the constriction, (b), (c), and 
(d): Distribution of C values obtained using purely electrical 
measurements, for (b) the same set of devices shown in Fig. 
[T] and (c),(d) for 310 additional devices on a separate chip, 
with active areas of both (c) 3 x 3.3 \im and (d) 10 x 10 \im. 
These data were obtained from 7g bs _R n , with their normaliza- 
tion set by single measurements of Lk(^bias) like that shown 
in (a), one for each of the three sets of devices shown in (b), 

(c) , and (d). The devices from (b) and (c) are nominally iden- 
tical in all respects, though they were fabricated on different 
NbN films, and at different times. Their C distributions are 
are therefore quite similar. The difference between (c) and 

(d) indicates that larger devices have a significantly higher 
probability of a sizeable constriction. 

Figs. MJd) and (c). Since optical testing of large num- 
bers of detectors is significantly more difficult and time- 
consuming than electrical testing, this is potentially an 
important screening technique. An example of data ob- 
tained with this technique is shown in Figs. [DJb), (c), 



and (d). For three sets of devices, we obtained C values 
from 7£ bs i? n , normalized using a single absolute measure- 
ment of C for each set via ik(^bias) as described above. 
Panel (b) shows data for the set of devices from Fig. [3a); 
panels (c) and (d) show results for a sample of 310 ad- 
ditional devices on a single additional chip, all having 90 
nm wide wires on a 200 nm pitch, but with active areas 
of both (c) 3 x 3.3 u.m, and (d) 10 x 10 |J.m. The devices 
from (b) and (c), which are nominally identical, exhibit 
very similar C distributions (though some yield fluctu- 
ation, which we commonly observe, is evident). Panels 
(c) and (d), however, clearly show a lower average C 
for the larger-area devices, qualitatively consistent with 
some fixed area density of constrictions. This implies 
that at the present state of film growth and nanowire 
patterning, the yield for high-DE single devices or device 
arrays covering 10 x 10 |J.m areas or larger will likely be 
quite low. 

As a final note, we remark on the obvious question 
of the origin of these constrictions. The most natural 
explanation would be lithographic patterning errors; for 
example, a small particle or defect in the resist before 
exposure could result in a localized narrow section of 
wire. However, we have performed extensive scanning 
electron microscopy of devices that were measured to be 
severely constricted (e.g. C ~ 0.5) and no such errors 
were observable. This suggests that constrictions in our 
devices result either from thickness variations or material 
defects, which may have been present in the film before 
patterning, and may even be due to defects present in 
the substrate itself before film growth. 

In conclusion, we have verified both optically and elec- 
trically that the large variations in detection efficiency 
between nominally identical superconducting nanowire 
single photon detectors are the result of localized con- 
strictions which limit the device current. Further work is 
ongoing to pin down the exact source of these constric- 
tions, with the hope of eventually fabricating large arrays 
of these detectors. 

This work is sponsored by the United States Air Force 
under Air Force Contract #FA8721-05-C-0002. Opin- 
ions, interpretations, recommendations and conclusions 
are those of the authors and are not necessarily endorsed 
by the United States Government. 
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